
use corrupt.dta,replace 
rename cons consume
tabstat  regist cost exp gdp ///
     indus inincome inout manage pop tax consume ///
	 ,stat(mean  sd)
	 

egen idd = group(id)
xtset idd year

xtreg inincome regist exp gdp indus inout pop tax consume,fe
est store m11
xtreg inincome regist exp gdp indus inout pop tax consume,re
est store m12
hausman m11 m12
*选择固定效应

xtreg cost regist exp gdp indus inout pop tax consume,fe
est store m21
xtreg cost regist exp gdp indus inout pop tax consume,re
est store m22
hausman m21 m22
*选择固定效应

xtreg manage regist exp gdp indus inout pop tax consume,fe
est store m31
xtreg manage regist exp gdp indus inout pop tax consume,re
est store m32
hausman m31 m32
*随机效应

esttab m11 m22 m31, compress replace scalar(r2 r2_a N F)  ///
      star(* 0.1 ** 0.05 *** 0.01) 

*取对数
gen lncost = ln(cost)
gen lnmanage=ln(manage)
global vars "cost manage inincome regist exp gdp indus inout pop tax consume"
foreach v of varlist $vars{
      gen ln`v' = ln(`v')
		}
	
	
xtreg lninincome lnregist lnexp lngdp lnindus lninout lnpop lntax lnconsume,fe
	est store m11
xtreg lninincome lnregist lnexp lngdp lnindus lninout lnpop lntax lnconsume,re
est store m12
hausman m11 m12
*选择固定效应

xtreg lncost lnregist lnexp lngdp lnindus lninout lnpop lntax lnconsume,fe
est store m21
xtreg lncost lnregist lnexp lngdp lnindus lninout lnpop lntax lnconsume,re
est store m22
hausman m21 m22
*选择固定效应

xtreg lnmanage lnregist lnexp lngdp lnindus lninout lnpop lntax lnconsume,fe
est store m31
xtreg lnmanage lnregist lnexp lngdp lnindus lninout lnpop lntax lnconsume,re
est store m32
hausman m31 m32
*固定效应

esttab m11 m21 m31 m32, compress replace scalar(r2 r2_a N F)  ///
      star(* 0.1 ** 0.05 *** 0.01) 
